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ABSTRACT 


We  study  a  family2of  steady  flows  of  an  ideal  fluid  in  a 
bounded  domain  D  c  R^  using  a  variational  method.  Adapting  an 
approach  due  to  Arnold,  we  characterize  dynamically  possible 
flows  as  constrained  extremals  of  the  kinetic  energy  considered 
as  a  functional  of  the  vorticity  «(x),  xeD.  We  restrict  our 
attention  to  solutions  having  the  special  form: 

“  -  ■»  -  1  1{*>0)  '  lD“\ax  “  1 

where  X  is  a  free  parameter  and  #  (the  streamf unction)  is 
defined  by  -A#  ■  u  in  D,  #  *  -W  on  3D  (u  is  an  undetermined 
positive  constant).  A  simple  proof  of  existence  of  solutions  is 
provided  by  a  direct  variational  method.  The  main  results 
concern  the  singular  limit  X  ♦  •  (which  implies  y  ♦  •  ) .  We 
prove  that  ♦  6  *,  the  Dirac  delta  measure  at  x*  e  D,  in 
the  distributionalxsense,  and  we  characterize  the  point  x*  in 
terms  of  the  geometry  of  D.  Further,  we  show  that  supp  u>^ 
tends  asymptotically  (in  a  strong  sense)  to  an  inf initesimaidisc 
about  x*.  This  follows  since  appropriately  scaled  version  of 
the  solutions  tend  to  the  unique  radial  solution  of  a 

corresponding  limit  problem.  In  fluid  dynamical  terms,  the 
classical  point  vortex  is  obtained  as  the  limit  of  vortices  of 
finite  cross-section  (of  the  Rankine  type). 
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SIGNIFICANCE  AND  EXPLANATION 


The  purpose  of  this  paper  (and  its  sequel)  is  to  give  a 
mathematically  rigorous  -  as  well  as  physically  natural  - 
discussion  of  certain  steady  solutions  of  the  Euler  dynamical 
equations  for  an  ideal,  two-dimensional  fluid.  The  flows 
considered  have  a  prescribed  distribution  of  vorticity 
«  =  curl  u  (u  devotes  thb  velocity  field),  and  are  separated 
into  regions  where  *  0  arid ■— 4»)  >  0.  The  shape  and  position 
of  the  'vortex  core*  (region  where, /£)  >  0)  for  a  flow  satisfying 
the  dynamical  requirements  is  then  determined  by  the  geometry  of 
the  fluid  domain  (assumed  bounded  here).  Solutions  of  the  fluid 
dynamical  equations  are  most  conveniently  characterized  by  a 
variational  principle  which  involves  finding  an  extreme  value  for 
the  kinetic  energy  of  the  flow  subject  to  certain  natural 
constraints.  This  approach  permits  a  precise  analysis  of  the 
properties  of  solutions  to  be  carried  out  in  a  unified  manner. 

In  this  respect  , special  emphasis  is  placed  upon  deriving  the 
(classical)  point  vortex  as  the  limit  of  solutions  with 
concentrated  vorticity — - 

The  occurence  of  vortices  (or  "eddies")  in  real  fluid  flows 
is,  of  course,  well  known.  The  traditional  literature  is  most 
often  restricted  to  idealizations  in  which  (for  two  dimensions) 
the  vorticity  is  concentrated  (as  Dirac  deltas)  at  a  finite 
number  of  points.  Also,  numerical  simulations  of  vortex  flows 
are  frequently  based  upon  this  discrete  vortex  approximation. 

The  present  work,  however,  is  devoted  to  a  simple  model  problem 
for  which,  first,  the  appropriate  mathematical  theory  is 
developed  without  the  above  mentioned  restriction,  and  second, 
the  (asymptotic)  nature  of  the  traditional  idealization  is 
analysed  qualitatively.  It  should  also  be  remarked  that  some  of 
the  methods  developed  here  for  steady  flows  can  be  adapted  to  the 
corresponding  time-dependent  flows. 
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existence  of  solutions  and  to  determine  the  qualitative 
properties  of  these  solutions.  We  first  give  a  simple 
proof  of  existence  using  a  direct  (variational)  method, 
and  compute  the  somewhat  novel  variational  conditions 
satisfied  by  extremals.  We  then  turn  to  our  principal 
results  which  concern  the  asymptotic  analysis  of 
solutions  (each  representing  a  vortex  of  finite  cross- 
section)  in  a  certain  singular  limit  (the  limit  being  a 
classical  point  vortex).  In  the  latter  analysis  we 
feel  that  an  advantage  is  achieved  here  over  other 
approaches  to  this  and  related  problems  since  sharp 
results  on  the  asymptotic  nature  of  solutions  are 
obtained  by  relatively  elementary  means. 

We  confine  ourselves  in  the  present  paper  (Part  I) 
to  a  very  simple  flow  geometry  :  the  fluid  domain  is 
bounded  and  the  flow  is  everywhere  tangential  on  the 
boundary.  In  a  subsequent  paper  (Part  II)  we  shall 
apply  the  results  developed  here  in  the  simple  case  to 
other  more  complicated  situations.  In  particular,  we 
shall  consider  steady  vortex  pairs  and,  especially, 
steady  vortical  wakes  behind  a  symmetric  obstacle  in  a 
uniform  stream.  Axisymmetric  analogues  (such  a  vortex 
rings)  can  also  be  studied  using  similar  methods. 

41.  VARIATIONAL  PROBLEM 

2 

Let  D  c  R  be  a  bounded  and  simply-connected 


domain  with  a  smooth  boundary,  3D.  Throughout  the 

sequel  we  shall  assume  for  convenience  that  diam  D  *  1; 

this  is  achieved  after  appropriately  scaling  the 

spatial  variables,  which  we  shall  write  as 

2 

x  =  (x^Xj)  e  R  . 

We  consider  a  steady  flow  of  an  ideal  fluid  of 
unit  density  in  the  (fluid)  domain  D.  The  velocity 
field  u(x)  *  (u^x),  u2(x))  and  the  pressure  p(x) 
are  required  to  satisfy  the  Euler  equations: 

(1.1)  7»u  =  0  in  D  and  v»u  ■  0  on  3D 

(1.2)  u*7u=  -7p  in  D  , 

where  v  is  the  exterior  unit  normal  on  3D.  The 
kinematical  conditions  (1.1)  are  satisfied  whenever 
there  exists  a  stream  function  i|»*i|>(x),  xeD,  in 
terms  of  which  the  velocity  field  can  be  represented  as 

(1.3)  u  =  J7i|i  ■  (♦  ,-i|>  )  in  D  and 

x2  X1 

iji  =  const,  on  3D? 

here  and  in  the  sequel,  J(a^,a2)  =  (-aj.a^)  denotes 
clockwise  rotation  through  n/2.  We  now  express  the 
dynamical  condition  (1.2)  in  a  weak  form,  which  will, 
in  turn,  lead  to  its  variational  characterization.  The 
vorticity  w  *  u>(x),  x  €  D,  is  defined  to  be 

(1.4)  a)  =  u2,x  -  u1#x  =  -Ai|>  in  D  . 

It  follows  by  a  well-known  calculation  that  then  (1.2) 


becomes 

(1.5) 


o)7ij)  =  -7  (p+  ^2  I  7i|»  |  ^  )  in  D  . 


We  recall  that  a  vector  field  a  =  a(x)  can  be 


expressed  as  a  gradient  a  =*  Vb  for  some  function 

b  =  b(x)  if  and  only  if 

/  a*JV$dx  =  0  for  all  t  e  c!!(D)  . 

D  u 

Applying  this  to  (1.5)  we  have  that  a  pressure  p 
exists  if  and  only  if 

(1.6)  /  <od  (♦,♦  )dx  =  0  for  all  +  e  C~(D)  ; 

D  u 

here  we  use  the  shorthand  notation 


*(*,♦)  -  -  *  A  -  *  ♦, 


We  let  the  Green  function  for  -A  in  D  with 
homogeneous  Dirichlet  data  on  3D  be  written 

(1.7) g(x,xa)  -  igi  loglx-x'l-1  -  h(x,x‘)  x,x'  e  D  . 

The  corresponding  Green  operator  we  denote  by 

(1.8)  Ge(x)  -  /  g(x,x,)»(x’)dx'  . 

D 

We  also  define  (for  later  use)  the  Kirchhof f-Routh 
function  to  be 

(1.9)  H(x)  »^h(x,x)  x  e  D  . 

In  terms  of  these,  equation  (1.4)  inverts  to  become 

(1.10)  ♦  -  G»  -  y  #  **  ■  -  ♦Ijp  . 

Returning  to  (1.6)  the  dynamical  requirement  is  now 
expressed  as  a  condition  on  u  alone: 

(1.11) /  u>3(Gu>,$)  dx  *  0  for  all  $  e  c!!  (D)  . 

D  u 

Thus,  the  fluid  dynamical  problem  (1.1, 1.2)  is  reduced 

to  the  problem:  find  w  such  that  (1.11)  holds.  We 
proceed  by  characterizing  such  u>  as  the  constrained 


extremals  of  a  certain  energy  functional. 

The  kinetic  energy  of  the  flow  is 

E  =  V2  /  I  u  |  2dx  =  V2  /  1 7(3W  I  2<*x  f 
D  D 

so  E  defines  a  functional  of  w,  which  upon 

integration  by  parts  becomes 

E (<o  )  =  V2  /  “(x)Go»(x)dx 
D 

(1.12) 

«  V2  /  /  g(x#x*  )a»(x)u»(x' )dxdx*  . 

D  D 

We  now  construct  a  family  of  variations,  of  w 

for  any  given  ♦  as  in  (1.11).  Let  x  *  $t(y)  denote 
the  solution  of  the  initial  value  problem 
||  *  JV*(x)  ,  x(0)  —  y  6  D 

in  a  small  interval  t  e  ( -T , T ) ,  say.  The 
transformations  y  ♦  5t(y)  then  form  a  (smooth)  one- 
parameter  family  of  area-preserving  dif feomorphisms 
of  D.  The  variations  are  defined  by 

o>(t)(x)  =  «(5”1(x))  x  e  D  . 

A  calculation  employing  the  symmetry  property 
g(x,x')  =  g(x',x)  then  yields  the  first  variations 

E(m^)  »  V2  /  /  g(x,x'  ^U”1^)  )<oU^1(x' )  )dxdx' 

D  D 

=  V2  /  /  gUt(y),  5t(y'  )  )u)(y)aj(y'  )dydy' 

D  D 

=  E(w)  +  t  /  u)3  (Gu>,  $  )dx  +  o(t) 

D 


-5- 


as  t  ♦  0.  Consequently,  (1.11)  is  equivalent  to 

(i.i3)  0  =  it  I t=o  * 

Therefore,  we  have  shown  that  (local)  extremals  of 
E(<i>)  over  any  class  of  admissible  functions  w 
chosen  wide  enough  to  contain  all  variations  of  the 
form  will  yield  solutions  of  the  original 

problem  (1.1,  1.2)  -  that  is,  dynamically  possible 
steady  vortex  flows  in  D. 

A  natural  choice  for  the  class  of  admissible 
functions  for  E(u)  is  the  class  of  (measure- 
theoretic)  rearrangements  of  a  given  function  u>0; 
that  is,  all  those  u  satisfying 

meas  {x  €  D  :  u(x)  >  X} 

*  meas  {xeD  :  o>q(x)  >  X}  for  all  X  e  R  . 

At  least  formally,  any  (local)  extremal  for  E(w)  over 
this  class  of  rearrangements  yields  a  solution  u  of 
(1.11)  having  the  prescribed  measure  distribution 
belonging  to  Wq.  A  variational  principle  of  this  form 
has  been  proposed  by  Arnold  Cl]  (Appendix  2)  and 
certain  criteria  for  the  stability  of  steady  flows  have 
been  derived.  Also  Benjamin  [2]  has  discussed  an 
existence  and  stability  theory  for  steady  vortex  rings 
in  this  context.  In  the  present  work  we  choose  to 
consider  instead  a  wider  class  of  admissible  functions 
in  order  to  both  facilitate  the  (rigorous)  analysis  and 
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obtain  solutions  having  an  especially  simple  form.  We 
wish  to  model  a  steady  flow  in  D  which  is  separated 
into  an  irrotational  region  and  a  region  of  constant 
(positive)  vorticity;  we  do  not  permit  vortex  sheets 
interior  to  D.  To  this  end  we  let  K^(D),  the  class 
of  admissible  functions/  consist  of  all  those 
0  <  o>  e  l"(D)  satisfying  the  constraints 

(1.14)  /  a»(x)dx  =  1 

D 

(1.15)  ess  sup  w(x)  <  X 

x  e  D 

for  a  parameter  X  >  (meas  D)-1.  The  circulation  in 
(1.14)  is  normalized  to  be  1  after  scaling  the 
dependent  variable  to 7  X  plays  the  role  of  a  (vortex) 
strength  parameter. 

The  existence  result  proved  in  §2  is  the 
following:  there  exists  an  absolute  maximizer, 
o>  =  e  K^(D),  for  E  over  K^(D)  and  it  has  the 
special  form 

(1.16)  u>  *  xin,  =  {x  e  D  :  <|»(x)  >  0} 

where  ♦  is  defined  by  (1.10)  and  the  constant  u  >  0 
is  uniquely  determined  by  u»;  denotes  the 

characteristic  function  of  the  set  A.  The  statement 

(1.16)  represents  the  strong  form  of  the  variational 
conditions  for  the  extremal  u>  (it  implies  the  weak 
form  (1.6)),  and  y  arises  as  a  Lagrange  multiplier 


for  the  constraint  (1.14).  The  streamfunction  then 
solves  a  nonlinear  free  boundary  problem  which  can  be 
expressed  compactly  as  an  elliptic  equation  with 
discontinuous  nonlinearity* 

(1.17)  -A*  =  Xi{<)>0}  in  Dr 

we  note  that  the  velocity  field  u  =  JVi|»  is  continuous 
across  the  free  boundary,  3Q.  In  other  words,  when 
<»  takes  the  special  form  (1.16)  the  general  condition 
(1.6)  reduces  simply  to  the  requirement  that  dft  be  a 
streamline  ( for  ip ) . 

There  is  no  corresponding  uniqueness  result  for 
solutions,  or  even  for  absolute  maximizers  of  the  form 
(1.16).  Indeed,  for  certain  geometries  nonuniqueness 
can  be  proved  using  the  asymptotic  results  in  §3-4. 

The  main  body  of  analysis  occurs  in  §3-4  where 
the  asymptotic  nature  of  maximizers  as  X  ♦  •  is 

determined.  We  expect  to  obtain  in  the  limit  a  steady 
point  vortex;  indeed  we  show  that 

(1.18)  o>^(x)  *  fi(x-X*)  as  X  ♦  • 

where  the  convergence  is  in  the  sense  of  distributions 
and  6(z)  denotes  the  unit  (Dirac)  delta  measure 
concentrated  at  z  *  0.  Furthermore,  the  location 
X*  6  D  is  determined  (by  the  geometry  of  D)  according 
to 

(1.19)  H(X*)  =  min  H(x)  (recall  (1.9))  . 

x  e  d 


Strictly  speaking,  as  X*  is  not  necessarily  uniquely 
determined  the  statement  (1.18)  may  hold  only  for  a 


sequence  X  =  X^  In  §3  the  basic  asymptotic 

estimate  for  large  X  -from  which  all  else  follows  - 
is  proved;  it  states  that 

(1.20)  diam  (supp  u>^)  <  Re  (Xxe2=l) 

for  a  positive  constant  R  >  1  (independent  of  X). 

In  $4  we  establish  the  sharp  limiting  forms  of 

and  30^.  To  do  so  we  introduce  the  center  of 
vorticity 

X.  ■  /  x  (x)dx  , 

A  D 


and  the  scaled  version  or  ("blow-up")  of 

CA(y)  =  X-1ux(Xx  +  ey)  (Xife2^) 

»  iA  (y)  y  e  br(o), 

X 

where  A^  c  BR(0)  is  the  scaled  version  of  We 

prove  that  a  unique  radial  limit  is  approached  as 


B1(0) 


weakly  star  in  L  (D) 


3A^  +  313^0)  C1  sense  as  curves  . 

These  results  depend  upon  some  symmetry  lemmas  in 
potential  theory  given  at  the  beginning  of  §4. 

We  remark  that  the  condition  (1.19)  determining 
X*  is  anticipated  by  the  classical  theory  of  a  point 
vortex  moving  through  an  ideal  fluid  (otherwise 
irrotational)  within  a  bounded  domain.  According  to 
the  theory  developed  by  Ktrchhoff  and  Routh,  a  point 


vortex  of  unit  circulation  concentrated  at  X(t)  @  D 
(t  is  time)  moves  along  a  trajectory  determined  by  the 
system 

(1.21)  H  =  -JVH(X)  , 

evidently  a  Hamiltonian  system  with  one  degree  of 
freedom.  Now  we  see  why  the  limit  as  X  ♦  •  of  steady 
flows  induced  by  u>^  should  tend  to  a  steady  point 
vortex  at  X*  satisfying  (1.19);  indeed,  X*  is  an 
equilibrium  point  of  the  Hamiltonian  system  (1.21) 
since  VH(X*)  =  0.  For  a  full  discussion  of  the 
classical  theory  we  refer  the  reader  to  Lin  [12];  see 
also  the  recent  work  of  Richardson  [15]. 

It  is  possible  to  study  a  closely  related  class  of 
variational  problems  in  which  solutions  of  the  form 
(1.16)  are  replaced  by  solutions  of  the  following  forms 
u=X(4»  +  )®  0  <  8  <  "  (♦  as  above)  . 

The  appropriate  class  of  admissible  functions  is 
Kw(D)  c  L^D)  which  is  obtained  by  removing  the 
constraint  (1.15),  and  the  appropriate  functional  is 

Eg  (oi )  *  E (<u )  -  /  £  [“i*)]Pdx 
D 

for  p  *  1  +  1/8;  X  >  0  is  prescribed.  The  existence 
of  an  absolute  maximizer  u  =  u».  for  Ea  over 

A  #  P  -P 

KW(D)  can  be  proved  using  a  direct  method  similar  to 
that  given  in  §2;  general  existence  theorems  of  this 
type  appear  in  Berestycki  and  Brezis  [3],  Asymptotic 
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results  for  e.  as  X  *  »  can  also  be  shown  by 
a  ,  p 

natural  (though  not  entirely  obvious)  extensions  of  the 
techniques  developed  In  $3-4.  In  particular,  the 
estimate  (1.20)  for  large  \  holds  for  arbitrary 
0  <  0  <  •  (fixed);  we  will  not  present  the  proof 
here.  This  asymptotic  estimate  has  been  proved  already 
for  3  »  1  by  Caffarelli  and  Friedman  [6],  and  for 
0  >  1  by  Berger  and  Fraenkel  [53;  however  these 
authors  begin  with  a  different  variational 
characterization  of  solutions  from  the  above.  Their 
proofs  (fundamentally  the  same  in  both)  are  entirely 
different  from  our  proof  ($3).  First  the 
connectedness  of  0^  ^  ■  (♦  >  0}  must  be  established 
(this  seems  to  require  0  >  1  unless  special  symmetry 
of  the  domain  D  is  assumed),  and  then  the  (potential- 
theoretic)  capacity  of  g  in  D  is  estimated;  this 
yields  (1.20)  since  the  capacity  can  be  related  to  the 
diameter  of  g •  The  techniques  given  in  §3  seem 
preferable,  however,  because  (i)  they  are  simpler  and 
(ii)  they  are  more  closely  allied  with  the  techiques 
needed  to  study  analogous  time-dependent  flows  with 
concentrated  vorticity;  we  intend  to  pursue  this  latter 
topic  elsewhere.  Finally,  we  remark  that  the 
functional  E^  represents  a  penalization  of  the  energy 
functional  E  for  small  0,  and  in  fact  the 
solutions  w.  (as  in  (1.16))  can  be  obtained  as  the 


limits  of  the  (penalized)  solutions  o>x  ^  as 
0  ♦  0+.  Accordingly,  alternate  proofs  of  all  the 
theorems  of  §2  can  be  based  on  this  procedure. 

Related  papers  not  already  cited  include  - 
Goldstik  [9],  Keady  and  Norbury  [11],  Keady  [10], 
GallouSt  [8];  in  the  context  of  vortex  pairs  and/or 
rings  -  Berger  and  Fraenkel  [4],  Norbury  [14],  Friedman 
and  Turkington  [7].  The  applied  literature  is 
referenced  in  Saffman  and  Baker  [16]. 

§2. EXISTENCE 

We  recall  from  §1  the  definitions: 

KX(D)  =  {u»  e  L*(D)  :  /  u(x)dx  =  1  , 

(2.1)  D 

0  <  <a(x)  <  X  a.e.  x  e  D> 

(2.2)  E (<■»)  *  V2  /  /  g(x,x*  )»(x)w(x '  )dxdx'  . 

D  D 

We  assume  throughout  the  sequel  that 

(2.3)  X  >  (meas  D)”1  . 

An  absolute  maximum  for  E  over  KX(D)  is  now 
found  by  a  simple  direct  method. 

Theorem  2.1.  There  exists  w  =  6  K^(D)  such  that 

(2.4)  E(w)  =  ^  max  E(u>)  . 

u  e  kx(d) 

Proof.  We  shaLl  often  use  the  obvious  estimate  (recall 
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diam  D  *  1) 

(2.5)  g(x,x')  <  log  |x  -  x'T1  x,x'  e  D  . 

That  E  is  bounded  above  on  KX(D)  is  evident  from: 

2E(w)  <  sup  Gw(x) 
x  e  D 

<  sup  /  log  lx  -  x' |“1w(x' )dx' 
x  e  D  z  D 

<  4—  /  log  lyl-1dy  (X*e2=l) 

2*  lyl  «  e 

“  (log  ^-  +  ^  )  *  C(X)  . 

Thus  there  is  a  sequence  6  K^(D)  such  that 

lim  E(w  .)  »  sup  E(»)  <  •  . 

3  «  e  kx(d) 

4D 

As  KX(D)  is  clearly  a  compact  subset  o£  L  (D)  in 
the  weak  star  topology  we  can  extract  a  subsequence  of 
(call  it  again  w^)  such  that 

OD 

a>j  ♦  <i)  e  Kx(D)  weakly  star  in  L  (D)  . 

Now  it  remains  to  prove  E(«j)  ♦  E(w);  but  this 
follows  easily  since  g(x,x')  e  L1(DXD)  and 
“j(x)“j(x')  ♦  w(x)w(x’)  weakly  star  in  L~(DXD). 

This  completes  the  proof. 

The  calculations  of  §1  imply  that  any  maximizer 
u  yields  a  dynamically  possible  steady  flow?  we  state 
this  fact  next.  We  note  that,  by  standard  potential 
theory,  Gw  6  H2  (D)nH?:  (D)  and  Gw  e  for  every 


1. 


Corollary  2.2.  Whenever  <n  6  satisfies  (2.4)  then 

there  holds 

(2.6)  /  a>3(Gu>,<fr)dx  =  0  for  all  ♦  e  cj  (D)  . 


Integrating  by  parts  formally  in  (2.6)  we  get 
3(Ga>,u)  =  0.  This  means  that  there  is  a  functional 
dependence  F  (VF^O)  such  that  F(Gw,u>)  =  0 
identically.  Condition  (2.6)  may  be  said  to  express 
the  functional  dependence  between  Gw  and  u>  "weakly 
therefore.  We  now  give  this  dependence  explicitly. 


Corollary  2.3.  Whenever  w  e  (D)  satisfies  (2.4) 
then 

(2.7)  w  *  Xlft  a.e.  in  D,  fl  =  {x  e  D:  Gu>(x)  >  y} 
for  a  constant  y  >  0  uniquely  determined  by  o>;  lQ 
denotes  the  characteristic  function  of  0. 


Proof.  We  consider  a  family  of  variations  of  w 
different  from  that  used  in  the  derivation  of  (2.6); 
namely,  we  define 

w(s)(x)  *  w(x)  +  s  Czo(x)  -  Zl(x)]  s  >  0 

00 

for  arbitrary  z^.z^  e  L  (D)  satisfying 


/  zn ( x ) dx  -  /  z. (x)dx 
D  D  1 


Z0’Z1  >  0 

a.e.  in 

D 

zo  “  0 

a.e.  in 

D\(w  <  A  -  d) 

z^  ■  0 

a.e.  in 

D\ (w  >  6} 

for  some  5  >0.  When 

8  >  0  is 

sufficiently  small 

(depending  on  6,  lzQl- 

we  have 

w , a .  6  K^(D),  and  hence  we  conclude  that 

0  >  h  E(“(S))  ,B-0+ 

-  /  zn(x)Gw(x)dx  -  /  z. (x)Gu(x)dx  . 

D  u  D  A 

It  is  easy  to  see  that  the  latter  inequality  holds  for 
arbitrary  zg'zl  as  above  and  arbitrary  6  >  0  only 
if 

(2.8)  ess  sup  Got (x)  <  ess  inf  Gw(x)  . 

w (x) <A  w(x)>0 

The  continuity  of  Gw  implies  easily  that  strict 
inequality  cannot  hold  in  (2.8);  hence  we  may  define 

(2.9)  p  ■  ess  sup  Gw(x)  *  ess  inf  Gw(x)  . 

<i)  ( x )  <A  w(x)>0 

From  this  it  is  now  clear  that  w  =  A  a.e.  in 
{Gw  >  p)  and  w  ■  0  a.e.  in  {Gw  <  u).  Therefore, 
in  order  to  conclude  the  desired  representation  (2.7) 
it  suffices  to  show  that  w  *0  a.e.  inS  *  (Gw  *  p) 
(as  we  have  not  excluded  meas  S  >  0).  To  prove  this 
we  observe  that  both  Gw  and  VGw  when  considered  as 
functions  of  either  x^  or  x2  alone  are  absolutely 
continuous  on  almost  all  lines  (segments  in  D) 
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parallel  to  the  coordinate  axes.  Thus,  as  almost  all 
points  of  S  are  points  of  density  of  S  along  such 
lines,  it  follows  that  the  pointwise  first  and  second 
partial  derivatives  of  Gw  vanish  a.e.  in  S.  Hence 

w  =  -A Gw  *  0  a.e.  in  S,  as  required.  This 

reasoning  is  taken  from  Morrey  [13]  Theorem  3.2.2(c). 
Finally  we  notice  that  u  >  0  since  otherwise  w  *  X 
a.e.  in  D,  contradicting  hypothesis  (2.3).  This 
completes  the  proof. 

We  leave  the  general  question  of  the  regularity  of 
the  free  boundary  3{^>0}C(4»*0}  unresolved 
here.  All  subsequent  results  of  this  paper  depend  upon 

(2.7)  only.  In  §4,  however,  we  do  prove  that  the  free 

boundary  is  a  simple  closed  C1  curve  if  X  is  large 
enough.  For  arbitrary  X  we  conjecture  that  the  set 

of  irregular  free  boundary  points  (where  Vf  *  0)  is 

finite  or  has  one-dimensional  Hausdorff  measure  zero, 
but  this  appears  to  be  an  open  problem. 

§3.  ASYMPTOTIC  ESTIMATE 

This  section  is  devoted  to  the  proof  of  estimate 
(1.20)  which  is  basic  to  the  asymptotic  analysis  of 
solutions  as  X  -*•  « .  We  suppose  throughout  that  w 
satisfies  (2.4),  and  we  write  <J>  =  Gu>  -  y  according  to 

(2.7) .  For  any  such  w  we  define 
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(3.1)  T(«* )  ■  ^2  I  IV<|>  +  |2dx  (♦+  *  max  (♦ ,  0}  )  ; 

D 

T  is  the  kinetic  energy  of  the  vortex  core 

a  -  (♦  >  0} . 

We  prepare  the  proof  of  the  main  theorem  by  two 
short  lemmas.  We  write  C,  C^#C2«...«  for  positive 
constants  independent  of  X . 

Lemma  3.1.  There  holds  the  lower  bound: 

(3.2)  E(«)  >  ^  log  |  -  Cx  (Xwe2  -  1)  . 

Proof.  We  consider  for  large  X 

(3.3)  "  -  “b  «)  9  K»«D> 

A  ^  A 

where  X  e  D  is  chosen  such  that  H(X)  *  min  H(x) 

xeD 

(the  reader  will  verify  that  H(x)  +  •  as  x  ♦  3D). 

Then  we  may  estimate 
* 

E(u>)  >  E ( (o  ) 

■  T-  /  /  log|x-x'  |"1(i)(x)«i>(x' )dxdx'-H(X)  +o(l) 
D  D 

>  47  lo9  2~  -  H(X)  +  o(l) 
as  X  -*•  •.  From  this  (3.2)  clearly  follows. 

Lemma  3.2.  There  holds  the  upper  bound: 

(3.4)  T(w)  *  V2  /  ’f’wdx  <  C» 

D  1 


Proof .  Since  i|»  =  -y  <0  on  9D  we  may  integrate  by 
parts  to  get 


2T(w)  -  /  Vip  »V</dx  =  /  <|i  oadx  =  /  'I'udx, 

D  D  D 

using  (2.7)  in  the  last  equality.  Furthermore,  we  have 
/  i|>«dx  =  X  /  <>+dx  <  X  (meas  ft)  ^  {  /  (^+)2)  ^2  . 


Thus,  by  applying  the  Sobolev  inequality  to  4> 
{  /  (i|>+)  2dx)  <  c  /  |V*+|dx, 


we  then  obtain 


2T(»)  <  CA(meas  0)  ^  J  |Vi|i+|dx 


<  CX(meas  fl){  /  |Vi)»+|2  dx}  ^2 
D 


=  C  [2T(w) ]  y2  . 


this  establishes  (3.4). 


Theorem  3.3.  There  is  a  constant  R  >  1 
( independent  of  X )  such  that 

2 

(3.5)  diam  (supp  u> )  <  Re  (Xire  =1)  . 


Proof .  We  estimate  the  parameter  m  using  the  obvious 


identity 

(3.6) 


T(w  )=  e(w  )  -  V2  u  ; 


then  Lemmas  3.1  and  3.2  combine  to  give 


3. 


2 


Now  let  x  6  supp  o>  be  fixed  arbitrarily.  Then 
Gw (x)  >  v ,  and  so  we  find  (recalling  (2.5)) 

57  1o9  7  -  C3 

<  /  g(x,x*  )o)(x'  )dx*  <  *7  J  log  |  x-x  '  |_1o) (x '  )dx ' 

D  D 

Equivalently,  we  write  (for  any  R  >  1) 

-  2wC3  <  /  log  «(x'  )dx* 

-  /  +  /  log  i  e  ■  |  o)  ( x ' )  dx '  . 

B__  (x)  D\B„  (x)  |X  X  1 


We  observe  now  that 

/  log  i/.tt  «(x'  )dx'  <1  /  log  t~t  dy  =  ^  . 

BRe(x)  lyl<e  7 

Therefore  we  get 

-2*C--16<  /  log  .  l  ,  ■  ***  ( x 1 )  dx 

3  D\B„(x)  |X  X  1 


‘  (1°9  5>  DVB^(x)“U')ax‘  ' 


after  manipulation  this  becomes 


(3.8) 


/  u>(x')dx'  <  C.(log  R) 

d\BR£(x) 


We  now  see  that  the  desired  estimate  (3.5)  follows 

immediately  from  inequality  (3.8).  Indeed,  we  claim 

that  diam  (supp  m)  <  2Re  if  R  is  fixed  large  enough 

to  ensure  that  C^(log  R)  ^  <  V2  »  Otherwise  there 
1  2 

would  exist  x  ,  x  6  supp  u>  with  the  property  that 
B^dx1)  n  B_  (x2)  =  and  so,  by  (3.8), 


/  o)(x')dx'  >  /  ,  bt(x')dx'  +  /  0  w(x')dx'  >  1  ? 

D  BRe<x  >  BRe<x  > 

this  contradicts  the  constraints  for  u)  6  K^(D),  and 

so  the  theorem  is  proved. 

Remark .  Of  course,  the  constant  R  depends  upon  the 

domain;  an  examination  of  the  proof  shows  that  R  can 

be  taken  in  the  form:  R  *  c  exp  (kH* )  where 

H*  =  min  H  and  c,k  are  absolute  positive  constants 
D 

§4. LIMITING  BEHAVIOR 

In  this  section  we  study  the  limiting  form  of  the 
maximizers  w  =  u>^  supplied  by  Theorem  2.1  as 
A  ♦  ®.  We  find  that  their  appropriately  scaled 
versions  tend  to  a  unique  limiting  function  which  is 
necessarily  radially  symmetric.  To  this  end  we  first 
prove  some  symmetry  lemmas  in  potential  theory  which 
are  themselves  of  independent  interest. 

Lemma  4.1.  Suppose  that  an  open  set  A  with 

A  £  Br(0)  ,  0  <  R  <  ®,  possesses  the  property  that 

A  =  {x  6  Br{0)  j  V(x)  >  y)  for  some  constant 

Y  =  Y (A)  where  V  =  VA  is  the  potential  of  A: 

(4.1)  V(x)  =  /  log  | x-x ' | _1dx '  . 

^  A 

Then  A  is  necessarily  a  ball  (and  so  V  is  radial). 


Proof.  We  shall  utilize  a  reflection  (or  "folding") 
argument  very  similar  to  an  argument  used  by  Serrin 
[17]  for  a  different  symmetry  result  in  potential 
theory.  In  particular,  we  shall  show  that  A  is 
symmetric  across  a  line  x±  ■  const.;  this  suffices  to 
prove  the  radial  symmetry  of  A  since  we  may  rotate 
coordinate  axes  arbitrarily  without  changing  the 
conclusion. 

We  define  the  open  sets,  for  -R  <  t  <  R, 

At  *  A  n  {x  :  x1  <  t).  A*  ■  {x  s  (2t“X1#x2)  e  At)  ; 
that  is.  A*  is  the  reflection  of  At  across  the 
line  Xj_  ■  t.  We  then  consider  s  »  max  {t  :  A*  c  a) 
(it  is  easy  to  see  that  this  maximum  is  achieved  and 
that  -R  <  s  <  R) .  We  now  claim  that  A  is  symmetric 
across  the  line  Xy*  s;  we  intend  to  establish  the 
claim  by  showing  that 

(4.2)  meas  N  *  0  for  N  —  A\  (A  u  a*)  . 

S  8 

First  we  observe  that  there  must  exist  a  point 

x*  e  3A  n  3a*,  since  otherwise  A*  c  A  and  hence 
s  s  — 

A*  c  A  for  some  t  >  s  (contradicting  the  maximality 

of  s).  We  then  write  x  =  (2s-x?,xS)  e  3A  .  We  now 

i  z  a 

proceed  to  consider  the  two  distinct  cases: 

(1)  x^  <  s  <  x^,  (2)  x ^  =  a  =  x^.  The  reader  may 

check  by  example  that  either  case  may  occur. 

Case  1.  Since  x,  x*  e  3A  there  holds 


owing  to  the  fact  that  the  corresponding  integrals 
over  A_  and  A*  cancel  each  other.  But  the  above 

s  S 

integrand  is  strictly  positive  on  N\{x^  =  s}  and  so 
it  is  not  possible  that  meas  N  >  0;  thus  (4.2)  follows 
in  this  case. 

Case  2.  Whenever  x  =  x*  6  3 A  n  a a*  n  {x,=s}  there 

s  l 

holds 

vxx(x)  ■  hr  /  lx  ■  *'  r2(x‘-x1)dx*  , 

because  again  the  corresponding  integrals  over  As  and 

A*  cancel  each  other.  Thus  if  meas  N  >  0  then 
s 

V  (x)  >  0  (as  the  integrand  is  strictly  positive  on 
1 

nUx^s}).  Also,  using  the  fact  that  now  W  ^  0  at 
any  such  point  x  =  x*,  we  have  by  the  implicit 
function  theorem  that  3A  is  a  C1  curve  in  a 
neighborhood  of  x  =  x*.  But  then  it  is  not  difficult 
to  see  that  these  facts  together  imply  that  there  must 
exist  t  >  s  such  that  A*  C  a,  a  contradiction;  thus 
(4.2)  also  follows  in  this  case. 

Lemma  4.2.  Let  the  class  K*  consist  of  all  those 
0  <  c  e  L* (Br(0)),  1  <  R  <  «,  satisfying  the 

constraints 

o\  /  C  ( x)  dx  =  if  ,  ess  sup  c(x)  <  1  . 

'4•J,  Br(0)  x  6  BR( 0 ) 

Let  the  functional  F  be  defined  by 
(4.4)  F(? )  = 

4r  /  /  log  lx  -  x '  r1^  (x)c  (x’  )dxdx '  . 

4x  BD (0 )  Bo(0) 


1 

I 


Then 

over 

(4.5) 


c*  = 


IB1(0) 


is  the  unique  maximimizer  of 


K*  for  which 


/  x? (x)dx  =  0  . 

V°> 


p 


Proof.  The  methods  of  §2  apply  here  without 
change.  We  remark  that  the  existence  of  a  maximizer 
follows  by  the  proof  of  Theorem  2.1.  Furthermore#  the 
proof  of  Corollary  2.3  shows  that  any  maximizer 
C  of  F  over  K*  has  the  form 

C  *  IA  ,  A  -  {x  e  Br(0)  :  V(x)  >  Y) 
for  V  given  by  (4.1)  and  some  constant  y 
(determined  by  ?  as  before).  Therefore,  Lemma  4.1 
applies  and  we  conclude  that  A  is  necessarily  a  ball 
of  radius  1?  obviously,  (4.5)  is  imposed  to  determine 
the  center  of  the  ball  uniquely. 

Remark .  Both  of  the  preceding  lemmas  have  natural 
generalizations  to  n  dimensions  which  can  be  proved 
using  the  same  methods.  Also  a  wide  class  of  kernels 
(defining  the  potential  V  and  the  functional  F)  can 
be  permitted  in  place  of  the  Newtonian  kernel.  We 
leave  these  extensions  to  the  interested  reader. 

We  now  proceed  to  determine  the  limiting  behavior 
of  u>^  as  X  -*•  00 .  We  define  the  center  of  vorticity 
to  be 


(4.6)  X,  =  /  xcu.  (x)dx  . 

A  D 

We  now  fix  (for  the  remainder  of  the  discussion)  a 
sequence  X  =  X ^  ♦  »  such  that  there  exists  a  limiting 
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center 


(4.7)  Xx  ♦  X*  e  D  as  A  =  . 

Theorem  4.3.  Any  X*  as  in  (4.7)  satisfies 

(4.8)  H(X* )  =  min  H(x)  . 

X6D 

A  A 

Proof.  Let  w  =  a»x  be  defined  by  (3.3)  with 

A 

H(X)  =  min  H(x).  It  is  a  consequence  of  Lemma  4.2 
xeD 

that  then 

V2  /  /  h(x,  x'  )o»x  (x)wx  (x‘  )dxdx' 

D  D 

-  i—  f  f  log  |  x-x*  |  (x)w.  (x‘ )dxdx*  -  E ( u>  ) 

4ir  D  D  A  A  A 

<  i—  /  /  log  |x-x*  (x)a».  (x’  )dxdx’  -  E(o»  ) 

4*  D  D 

<  V2  I  I  h(x,x' )»x(x)«x(x' )dxdx'  , 

D  D 

A 

using  simply  E(u>x)  >  E(o>x)  last  inequality. 

Taking  A  ♦  •  in  the  above  we  obtain,  by  virtue  of  the 

A 

estimate  (3.5),  the  inequality  H(X*)  <  H(X);  this  is 
the  required  result  (4.8). 

We  can  now  say  that  (due  to  (3.5))  as  A  =  A  ^  « 

(4.9)  o>x(x)  ♦  6(x-X*)  as  distributions 

where  5(x-X*)  is  the  unit  (Dirac)  delta  measure  at 
x  =  X*  given  according  to  (4.8).  The  precise 
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asymptotic  nature  of  u»x  is  expressed  in  terms  of  its 
scaled  version  in  the  following  theorem. 


Theorem  4.4.  Let  Cx  e  L~(Br(0))  be  defined  by 

(4.10)  Cx(y)  *  X"1wx(Xx+ey)  (Xire2=l) 
for  fixed  R  as  in  (3.5).  Then  as  X  ♦  " 

(4.11)  cx  ♦  c*  -  IB  weakly  star  in  L~(BR(0))  . 


Proof. 

C  e  K*, 

<*>(x) 


We  retain  the  notation  of  Lemma  4.2.  For  any 
let  u  e  KX(D)  be  defined  as 


x  e  bR£(xx) 

X  e  D\BRe (Xx)  . 


A  direct  calculation  then  yields  as  X  • 
(4.12)  E(S)  =  ^  log  |  F(c) 


-  V2  I  J  h(x,x' )u(x)u(x' )dxdx* 
D  D 


~  log  i  “  H(X*)  +  o(1> 


and  likewise  this  statement  holds  for  w,  c  replaced 

by  u>x ,  ?x  .  But  then  E(«)  <  E(wx)  implies  that 

F(?)  <  F(CX)  +  od)  as  X  -►  ».  For  any  sequence 

Xn  >  •  there  is  a  subsequence  X^  ♦  »  such  that 

Cx’  ■►5  e  K*  weakly  star  in  L*(BR(0)).  Then 
m 

F(?)  =  lim  F(c.  ,)  >  F{?)  for  all  C  8  K*, 


and  so  we 


conclude  by  Lemma  4.2  that  c  «  C*  3  IR  (the 

construction  (4.10)  ensures  condition  (4.5)).  Since 
this  conclusion  is  independent  of  the  sequence  Xm 
taken  (by  the  uniqueness  of  c*)  we  find  that  (4.11) 
holds  whenever  X  ♦  •. 

The  conclusion  of  the  preceding  theorem  now 
implies  the  convergence  of  the  scaled  versions  of  the 
corresponding  streamf unctions. 

Theorem  4.5.  Let  v^  e  C^(BRl(0))  be  defined  by 

(4.13)  v*(y)  “  w*X^XA+ey^  (Xire2=  1) 

for  fixed  R  as  in  (3.5)  and  R  <  R'  <  •.  Then  as 


A  -*•  * 

(4.14) 

VA  +  v* 

in 

C1^,  (0) 

where 

v*  is  given  by 

Ji 

(l-lyl2)  0  < 

lyl  <  1 

(4.15) 

v*(y)  «  v*( |y| ) 

,  .  _i 

lyl  <  00 

|2 

log  |y|  1  < 

Remark.  The  limit  function  is  also  expressible  as  the 
potential 

(4.16)  v*(y)  “^7/2  lo9  ly~y ' l”1C*(y’ )dy'  . 

R  9  i  « 

We  note  that  v^ ,  v*  6  H  (BR,(0))  nc  '  (Br,(0))  for 
every  0  <  a  <  1,  and  -Av^  3  -Av*  =  c*. 

v*  the  (classical)  Rankine  streamfunction. 


We  call 


Proof.  We  define  Vx  6  (^(R^)  by 

vx(y)  *37/ 2  1o9  ly~y‘ r1cx(y')dy'  . 

R 

Since  supp  Cx  C  BR(0)  and  0  <  ?x<  1  a.e.  in  BR(0), 
we  may  apply  standard  potential  theoretic  estimates 
(see  Morrey  [13])  to  conclude  that 

|VVx(y)|  <  C 

|WX  (y1)-VVx  (y2)  I  <  C l yX-y2 1  log  (14 — ) 

1  2  ly  -y  I 

for  all  y,  y ,  y  6  BR,(0)?  C  is  a  absolute 
constant.  It  follows  that  both  families  {Vx>  and 
{ VVX }  are  equicontinuous  in  BR,(0).  Since  (4.11) 
implies  the  convergence  Vx  ♦  v*,  VVX+  w*  pointwise 
in  Br,(0),  we  find  then  that 

(4.17)  Vx  ♦  v*  in  C1(Br1(0))  as  X  *  -  . 

2 

Let  Dx  =  {y  6  R  :Xx+ey€D)  denote  the 

scaled  version  of  D.  Then  vx  is  defined  on  Dx  and 

vx  1 3 D  ■  *  “1,WX*  We  observe  that  by  the 

X 

identity  (3.6)  and  the  established  upper  and  lower 
estimates  for  E(u>x  )  (along  with  the  bound  for 
T(u»x))  there  holds 

ux  =  log  ^-+0(1)  as  X  ♦  •  . 

Thus,  since  d/e  <  |y|  <  1/e  for  all  y  63DX  if 
d  *  ^2  dist  (X*,3D)  <  dist  (XX#3D),  we  conclude  that 
vx(y)  =  V2  log  J y ! ~ 1  +  0(1)  =  Vx(y)  +  0(1)  as  X  ♦  • 
for  all  y  63  Dx.  Now  since  we  have 

Mvx-Vx>  *0  in  Dx  and  Ivx-Vxl  <  C  (a  constant 
independent  of  X),  we  apply  the  interior  gradient 
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estimate  for  harmonic  functions  to  get 


sup  |Vv. (y)  -  VV. (y) |  <  Cc  . 

Y6Bri (0) 

Therefore,  vx“va  ♦  c*  =  const,  in  C1(BR,(0));  now 
recalling  (4.17)  we  find  that  vx  ♦  v*  +  c*  in 
C1(B_,(0)).  Finally,  we  see  that  c*  =  0  since 

IV 

meas  (vx  >0}  =  meas  {v*  >  0}  *  *  by  virtue  of  the 
fact  that  Cx  *  I(vx>0}  e  K*‘ 

Corollary  4.6.  Let  Ax  *  {y  e  BR,(0)  :  vx(y)  >  0). 

Then  Ax  tends  asymptotically  to  the  unit  ball  as 
A  ♦  •  in  the  sense  that 

(4.18)  Ax  *  (y  *  (r  cos  8,r  sin  6):  0  <  r  <  ax(8)} 
where  ax  e  C1[0,2»]  is  a  periodic  function  and 

dax 

ax  ♦  1,  ♦  0  uniformly  in  [0,2w], 

Proof.  This  is  an  immediate  consequence  of  Theorem  4 . 5 
by  the  implicit  function  theorem. 

We  end  with  some  expansions. 

Corollary  4.7.  The  following  asymptotic  expansions 
hold  as  A  ♦  •  : 

(4.19)  E(u»x)  *  ~  (log  j  -  H (X* )  +  o(l) 
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(4.21) 


log  i  -  2H(X*)  +  o(l) 


Proof.  Derivation  of  (4.19)  i  Applying  (4.12)  to 
and  using  F(c^)  ♦  F(c*)  we  see  that  it  suffices  to 
calculate  the  value  of  F(c*);  this  is  given  by 


F(C*)  *  \  j  v*(y)dy  »  . 

*  n  /At 


B^O) 


Derivation  of  (4.20)  :  Changing  variables  we  have 


2T 


<v  -  4 .  |Sd*  *  4 .  L  .  ""*'**  -  4 


BR(0) 


*  B1(0) 


Derivation  of  (4.21)t  Using  (3.2)  this  is  immediate. 
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